Sub-Poissonian shot noise in non-degenerate diffusive conductors 
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A theory is presented for the universal reduction of shot noise by Coulomb repulsion, which was 
observed in computer simulations of a disordered non-degenerate electron gas by Gonzalez et al. 
[Phys. Rev. Lett. 80, 2901 (1998)]. The universality of the reduction below the uncorrelated value is 
explained as a feature of the high-voltage regime of space-charge limited conduction. The reduction 
factor depends on the dimensionality d of the density of states, being close but not quite equal to 
1/d in two and three dimensions. 

PACS numbers: 72.70.+m, 72.20.Ht, 73.50.Fq, 73.50.Td 



The motivation for this work comes from a remarkable 
recent Letter |[] by Gonzalez, Gonzalez, Mateos, Pardo, 
Reggiani, Bulashenko, and Rubf on the "universality of 
the 1/3 shot-noise suppression factor in non-degenerate 
diffusive conductors" . Shot noise is the time-dependent 
fluctuation in the electrical current caused by the dis- 
creteness of the charge. In the last few years there has 
been a breakthrough in the use of shot-noise measure- 
ments to study correlation effects in diffusive conductors 
|^| . In the absence of correlations between the electrons, 
the current fluctuations 5I(t) around the mean current I 
arc described by a Poisson process, with a spectral den- 
sity P at low frequencies equal to Pp isson = 2el. Corre- 
lations reduce the noise below the Poisson value. 

The Pauli exclusion principle is one source of correla- 
tions, Coulomb repulsion is another. In a degenerate elec- 
tron gas with elastic impurity scattering the reduction 
is a factor one-third ||J4|]. This reduction is due to the 
Pauli principle. The remarkable finding of Gonzalez et al. 
was that Coulomb repulsion in a three-dimensional non- 
degenerate electron gas also gives P/Pp isson = ■§■ They 
argued for a universal physical principle behind the one- 
third reduction of the shot noise from elastic scattering, 
regardless of whether the origin of the electron correla- 
tions is quantum mechanical (Pauli principle) or classical 
(Coulomb repulsion). 

The significance of Ref. jj] has been assessed criti- 
cally by Landauer |J. While in a degenerate electron 
gas the one-third reduction is independent of the num- 
ber d of spatial dimensions, Ref. M finds P/Pp isson = \ 
for d — 2 — spoiling the supposed universality of the 
reduction factor in non-degenerate conductors. Still, it 
remains a remarkable finding that the ratio P/Pp isson 
has no dependence on microscopic parameters (such as 
mean free path I or dielectric constant k) or global pa- 
rameters (such as temperature T, voltage V , or sample 
length L), as long as one stays in the high- voltage, diffu- 
sive regime (eV ^> kT and L 3> I). The findings of Ref. 
[pj are based on a computer simulation of the dynam- 
ics of an interacting electron gas. Here we develop an 
analytical theory that explains these numerical results. 



We use the same theoretical framework as Nagaev 
Q used for the degenerate electron gas, namely the 
Boltzmann-Langevin equation || in the diffusion approx- 
imation (valid for L 3> I). A difference with Ref. Q] 
is that the kinetic energy e = |mn 2 now appears as 
an independent variable. (In the degenerate case one 
may assume that all non-equilibrium electrons have ve- 
locity v equal to the Fermi velocity.) The charge density 
p(r, e,t) = p + Sp, current density j(r, e,t) = j + 6j, and 
electric field E(r, t) = E + 5~E fluctuate in time around 
their time-averaged values (indicated by an over line). In 
the low-frequency regime of interest we may neglect the 
time-derivative in the continuity equation, 
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j = 0. 



(1) 



Current and charge density are related by the drift- 
diffusion equation 



1 = -D% + 63, 
or oe 



(2) 



with a(e) = e 2 D(e)v(e) the conductivity, D(e) = DqE 
the diffusion constant, and v{e) = VQe d / 2 ^ 1 the density 
of states. (The coefficients Do and vq are e-independent, 
assuming an energy-independent effective mass and scat- 
tering rate.) The function f(r,e,t) = pjev = / + 5f is 
the occupation number of a quantum state. (In equilib- 
rium, the mean / is the Fermi-Dirac distribution func- 
tion.) The "Langevin current" S3(r,e,t) is a stochastic 
source of current fluctuations from elastic scattering ||. 
Its first two moments are S3 = and 



5Ji(r, e, 1)51^', e', f ) = 2cr(e)/(r, e)[l - /(r, e)] 

^5 l , J 5{v-Y , )5{e-e')5{t-t l ). (3) 

The definition of a non-degenerate electron gas is / <C 1, 
so that we may ignore the factor 1 — / in this correlator. 

We need one more equation to close the problem, 
namely the Poisson equation 



or 



de (p - p cq ) 



(4) 
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with p eq the compensating charge present in the semi- 
conductor in equilibrium (equal to the charge density of 
the carriers prior to the injection from the contacts). 

The geometry we are considering is a disordered semi- 
conductor of uniform cross-sectional area A sandwiched 
between metal contacts at x — and x = L. We de- 
note by r± the d — 1 dimensional vector of transverse 
coordinates. (Only d = 3 is physically relevant, but we 
consider arbitrary d for comparison with the computer 
simulations.) The mean values of p, j, and E are inde- 
pendent of rj_, but the fluctuations are not. We define 
the linear charge density p(x,t) — J dr± J de p(r,e,t), 
the electric field profile E(x, t) = A^ 1 J dr± E x (r, t), and 
the currents I{t) = J dr± J dej x (r,e,t) and 5J(x,t) = 
J dr± J de 5J x (r 7 e 7 t). The total current / is independent 
of x at low frequencies because of the continuity equation, 
but the Langevin current 6 J is not so restricted. In view 
of Eq. (0), the two currents are related by 



P = 2 



dr i 



deD^- + ppE + SJ, (5) 
ox 



with p — ^deDo the mobility. (To write the drift term 
in the form ppE we have made a partial integration over 
energy and linearized with respect to the fluctuations.) 

A non-fluctuating voltage V is applied between the two 
metal contacts, with the current source at x = and the 
current drain at x = L. (The charge e of the carriers is 
taken to be positive.) For high V the charge injected into 
the semiconductor by the current source will be much 
larger than the charge p cq present in equilibrium. We 
will neglect p eq altogether. For sufficiently high voltages, 
when all the surface charge at x = has been injected 
into the semiconductor, the system enters the regime of 
space-charged limited conduction, characterized by the 
boundary condition 



E(x) = at 



0. 



(6) 



The mean charge and field distributions in this regime 
were studied extensively in the past Q), but apparently 
the shot-noise problem was not. We argue that the uni- 
versality of the computer simulations jlj is a consequence 
of the homogeneity of the boundary condition (|6j). In- 
deed, if the boundary condition would have contained an 
external electric field, then the effect of Coulomb repul- 
sion on the shot noise would have depended on the rela- 
tive magnitude of the induced and external fields and 
hence on the value of k. No universal reduction fac- 
tor could have resulted. This scenario stands opposite 
to that in the degenerate case. There the reduction of 
shot noise occurs at low voltages, in the linear-response 
regime, when the induced electric field can be neglected 
relative to the external field 

The zero-frequency limit of the noise spectral density 
is given by 



dt5I(0)SI{t). 



(7) 



To compute P we need to relate the correlator of the 
total current SI to the correlator of the Langevin cur- 
rent SJ. At non-zero temperatures, SI contains also a 
contribution from the thermal fluctuations of the charge 
at the contacts (Johnson-Nyquist noise ||). This source 
of noise may be neglected relative to the shot noise for 
eV 3> kT, and we will do so to simplify the problem. 
The most questionable simplification that we will make 
is to neglect the diffusion term (cx dp/dx) relative to 
the drift term (cx E) in the drift-diffusion equation (j^). 
This approximation is customary in treatments of space- 
charge limited conduction 0, but is only rigorously jus- 
tified here in the formal limit d — > oo (when the ratio 
p/D -> oo). 

We are now ready to proceed to a solution of the cou- 
pled kinetic and Poisson equations. We consider first the 
mean values and then the fluctuations. Combination of 
Eq. (||) (without the term p cq ) and Eq. (||) (without the 
diffusion term) gives for the mean electric field 



E 



■ dE 



dx pnA 



E(x) 



2Ix 
pnA 



1/2 



(8) 



where we have used the boundary condition (g) . The x 
dependence of the electric field is the celebrated Mott- 
Gurney law [||. The corresponding charge density has 
an inverse square-root singularity at x — jnj. The 
corresponding voltage V = Edx oc y/I, so that the 
current increases quadratically with the voltage. These 
are well-known results for space-charge limited conduc- 
tion §. _ 

Linearization of Eqs. (|4|) and (|5|) around the mean val- 
ues gives for the fluctuations 

ox dx pnA 

^SE(x,t)=x-^ r^ 61 ^- 5 ^- (9) 



/„ {2ipKAy/ 2 ' 

A non-fluctuating voltage requires J Q L SE dx — 0, hence 

81{t) = 3 £ ^- (l - yfxft) SJ(x, t). (10) 

Combination of this relation between SI and SJ with Eqs. 
(^) and (]?]) gives an expression for the shot-noise power, 



P : 



36,4 r dx 



L Jo L 



(l - V^/L) 2 f dea(s)f(x,e). (11) 



To evaluate this expression we need to know the mean 
occupation number / out of equilibrium. For this pur- 
pose it is convenient to change variables from kinetic en- 
ergy e to total energy u = e + e(j){x), with 4>(x) the mean 



2 



electrical potential. Since E — —d<f>/dx, the derivative 
d/dx + eEd/de is equivalent to d/dx at constant u. The 
kinetic equations ([!]), (||) in the new variables x,u take 
the form 



9j 
dx 

The solution is 



f(x,u) 



0, j 



-a\u 



ecj>(x)} 



dx' 



df 
dx ' 



AR(u) J x a[u - e<j){x')\ 



(12) 



(13) 



where we have imposed the absorbing boundary condi- 
tion f{L,u) — at the current drain. (At high volt- 
ages the charge injected into the semiconductor by the 
current drain can be neglected.) The factor R{u) — 
A^ 1 J Q L dx/a[u — e(j)(x)] is the resistance of the semi- 
conductor. The mean current is related to R by el — 
J du /(0, u)/R{u). The argument u — e<f>{x) of a may 
be replaced by eV(x/L) 3 ^ 2 in the high-F limit. Then 
a oc x 3d / 4 and 



du a[u — e<fi(x)]f(x, u) 



el 
T 



dx 



AelL 



(3d -4)4 



'(?) 

(!) 



3ci/4 



3d/4' 



(14) 



Substitution into Eq. rtli]) yields our final result 



P = 



lUel 
3d -4 
24e7~ 



3d 2 -I- 22d + 64 



5 (d + 2)(3d + 4)(3d + 8) 



(15) 



The ratio P/Pp isson equals 0.341 and 0.514 for d = 3 
and d = 2, respectively, within error bars of the fractions 
| and i inferred by Gonzalez et al. from their computer 
simulations 0. The proximity of these numbers to d _1 
appears to be accidental. Indeed, for large d we find that 
-P/-Ppoisson - ► The large-d limit is a rigorous re- 

sult, while the finite-d values are not because we have 
neglected the diffusion term in Eq. ([|). 

In closing, we comment on the universality of the re- 
sults and on their experimental observability. Concerning 
the universality, the dimensionality dependence has al- 
ready been noted Q . For a given d there is no dependence 
on material parameters, however, the shot noise does de- 
pend on the model chosen for the energy-dependence of 
the elastic scattering rate. We have followed the com- 
puter simulations jjj] in assuming an e-independent scat- 
tering rate. In a model of short-range impurity scattering 
one would have instead a rate proportional to the den- 
sity of states. This would change the energy dependence 
of the diffusion constant from D cx e to D oc e 2 ~ d / 2 . 
The shot-noise power remains unaffected for d = 2, but 
for d = 3 one obtains Q P/Pp isson = 0.407 — some 



20% above the value for an e-independent scattering rate. 
Concerning the experimental observability, the main ob- 
stacle is the tendency of electron-phonon scattering to 
equilibrate the electron gas at the lattice temperature. 
Then instead of shot noise one would measure thermal 
noise (modified for a non-Ohmic conductor |ll]]), that 
is not sensitive to correlation effects. Experiments in 
degenerate systems have succeeded recently in observ- 
ing shot noise by reducing the sample dimensions to the 
mesoscopic scale The same approach may well be 

successful also in non-degenerate systems. 

I dedicate this paper to the memory of my father, 
the physicist J. J. M. Beenakker (1926-1998). Discus- 
sions with Oleg Bulashenko, Eugene Mishchenko, Hen- 
ning Schomerus, and Gilles Vissenberg are gratefully ac- 
knowledged. This work was supported by the Dutch Sci- 
ence Foundation NWO/FOM. 



[1] T. Gonzalez, C. Gonzalez, J. Mateos, D. Pardo, L. Reg- 
giani, O. M. Bulashenko, and J. M. Rubf, Phys. Rev. 
Lett. 80, 2901 (1998). 

[2] For a review, see M. J. M. de Jong and C. W. J. 
Beenakker, in: Mesoscopic Electron Transport, edited by 
L. L. Sohn, L. P. Kouwenhoven, and G. Schon, NATO 
ASI Series E345 (Kluwer, Dordrecht, 1997). 

[3] C. W. J. Beenakker and M. Biittiker, Phys. Rev. B 46, 
1889 (1992). 

[4] K. E. Nagaev, Phys. Lett. A 169, 103 (1992). 

[5] R. Landauer, Nature 392, 658 (1998). 

[6] Sh. Kogan, Electronic Noise and Fluctuations in Solids 
(Cambridge University, Cambridge, 1996). 

[7] M. A. Lampert and P. Mark, Current Injection in Solids 
(Academic, New York, 1970). 

[8] This is true in the zero- frequency limit. At non-zero fre- 
quencies the induced electric field is of importance in 
the linear-response regime as well, see: Y. Naveh, D. V. 
Averin, and K. K. Likharev, Phys. Rev. Lett. 79, 3482 
(1997); K. E. Nagaev, Phys. Rev. B 57, 4628 (1998). 

[9] N. F. Mott and R. W. Gurney, Electronic Processes in 
Ionic Crystals (Clarendon, Oxford, 1940). 
[10] A more accurate treatment including the diffusion term 
would cut off the singularity at a value set by the charge 
density in the contact. As a result, the electric field in the 
high-F regime would not vanish at x = but extrapolate 
to zero at a point x — ~x c inside the contact. Corrections 
to the shot-noise power of order x c /L are estimated to 
be less than 10% in the computer simulations of Ref. jl[ . 
[11] Thermal noise in a conductor with a non-linear I— V 
characteristic is given by P = 4kT(V/I){dI/dV) 2 . Since 
/ oc V 2 in the space-charge limited regime, one finds P = 
8kTdI /dV — two times larger than the thermal noise in 
an Ohmic conductor. This result can be derived from the 
theory presented in the text, by invoking the thermal 
equilibrium condition df/de = —J/kT. Eq. (O) then 



3 



takes the form P = 36/ifcTL" 2 / L dx (1 - y/x/L) 2 p(x), 

which can be evaluated using Eq. (^). 
[12] H. U. Schomerus, E. G. Mishchenko, and C. W. J. 

Beenakker (unpublished) . 
[13] F. Liefrink, J. I. Dijkhuis, M. J. M. de Jong, L. W. 

Molenkamp, and H. van Houten, Phys. Rev. B 49, 14066 



(1994); A. H. Steinbach, J. M. Martinis, and M. H. De- 
voret, Phys. Rev. Lett. 76, 3806 (1996); R. J. Schoelkopf, 
P. J. Burke, A. A. Kozhevnikov, D. E. Prober, and M. J. 
Rooks, Phys. Rev. Lett. 78, 3370 (1997); M. Henny, S. 
Oberholzer, C. Stru nk, and C. Schonenberger, preprint 
( |cond-mat/9808042j ). 



4 



